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We estimate the sensitivity of lensing observables to the parameters in the lens model 
(isothermal sphere/Navarro-Frenk- White profile) parameters and to cosmological parame- 
l- ters. We find that the observables are primarily dependent on the lens model parameters, 

' while the dependence on cosmological parameters is minor (especially so for the dark en- 

^0 ' ergy parameters) . We demonstrate the lens model degeneracy by deriving both the projected 

mass density profile and the circular velocity profile. We also identify a possible source of the 
problem of fitting the averaged mass profile of CL 0024+1654 with the Navarro- Frenk- White 
profile. 



1. Introduction. Gravitational lensing has been extensively studied as a useful 



\ cosmological tool to probe the high redshift universe. In particular, it has been 

' used to provide limits on cosmological parameters. For example, the statistics of 

^ ■ gravitational lensing of QSOs by intervening galaxies provide a powerful tool to 

set constraints on the cosmological constant. Lensed images of distant galaxies in 
cluster, so called arcs or rings, may provide a bound on the cosmological constant or 
• even on the equation of state of dark energy. This is mainly because the distance 

relation D^s/Ds, where D^s is the distance from the lens to the source and Ds is 
O ' that from the observer to the source, has a strong dependence on the cosmological 

+l3 . constant. 

I However, results obtained regarding gravitational lensing also depend on the lens 

model. Hence, in using it as a tool to obtain bounds on cosmological parameters, we 
- must be careful about the uncertainties inherent in the lens model itself. If derived 

r> I values of gravitational lensing observables are sensitive to the lens model, constraints 

' obtained on cosmological parameters may have a meaning whose implication will only 

become clear after resolving issues concerning the lens model. Only recently have 
these uncertainties begun to be systematically studied for the case of gravitational 
lens statistics. " In this letter, we elucidate the problem involving the model 
uncertainties by estimating the sensitivity of observables to the parameters in the 
lens model and cosmological parameters. We consider the giant arc system observed 
in CL 0024+1654 as an example. As an aside, we point out that the NFW profile 
fit may not be excluded by the current data. 

2. Lens model and observables. Let rj and ^ be the source position and the 
impact position in the lens plane, respectively. Define a length scale in the lens 
plane and a corresponding length scale rjQ = ^qDs/Di. Here Dl and Dg are the 
angular diameter distance to the lens and source, respectively. Then in terms of the 
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dimensionless quantities x = ^/^q and y = r]/rio, the lens equation is given by 

mix) 

y = x-^^, 0-1 

X 

where m{x) = 2 {U{x')/ X!cr)x'dx' , and Ucr = Ds / ^t^GDlDis- Here Dis is the 
angular diameter distance between the lens and the source. The determinant of the 
Jacobian A = dy/dy: of the mapping Eq.(O-l) is calculated as 



/V dx \ x 

The critical curves for axisymmctric lenses are those for which dct^ = 0. Circles 
satisfying m/x^ = 1 are called "tangential critical curves", while those satisfying 
d{m/x)/dx = 1 are called "radial curves". 

For an isothermal model with a core, the mass profile is p(r) = cr^/27rG(r^ + r^), 
with a being the one-dimensional velocity dispersion and Vc the core radius. The 
radius of the tangential critical curve, ^e, is then given by 

,filTei + ec = 47ra'^, (0-3) 

where 9c = Vc/Dl. The circular velocity at radius r is 

„V) = ^i^ = 2<,^fl-!ita.-^V (0-4) 

r \ r rcJ 

For the Navarro-Frenk- White (NFW) model, the mass profile is given by 
p(r) = Ps {r/rs)~^ {{r/rg) + Taking = fs-, the lens equation becomes ^-"^^ 

^PsTsOix) 

y = x--= , (0-5) 



cr 



where g{x) is defined by 



for x < 1 and by 



g{x) = \n^ + ^===tan-\^—. (0-7) 



2 

for a; > 1. The circular velocity is 

r V l + r/TsJ 

The truncated isothermal sphere (TIS) model is a particular solution of the 
Lane-Emden equation that results from the collapse and virialization of a top-hat 
density perturbation. The mass profile is fit well by 
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where {A,a'^ , B,b'^) = (21.38,9.08,19.81,14.62). Shapiro et al. fo^^^j ^j^^^ this 
fitting formula is accurate within 3% over < r/rc ^ 30 for both the Einstein de 
Sitter model and low-density models (J7m ^ 0.3). The circular velocity is 

v^(r) = AirGport [A- B tan~^ - H tan"^ - 

V X ax b 

where x = r/vc- 

At this stage, it may be illuminating to estimate the sensitivity of the observ- 
able {6e) to the parameters in lens model and the cosmological parameters. Using 
Eq.(0-3), for the case of an isothermal sphere with a core, the sensitivity is evaluated 
around Vc = 0, = 0.3, f^A = 0.7, w = —1 as 

-r^ + 2 7.5 X 10-2— iH. + 2.9 x lO'^ — , 0-11 

9e 9e 9c (t ilM w 

where w is the equation of state of dark energy {w = —1 for the cosmological con- 
stant), and the source rcdshift and the lens redshift are taken as Zg = 1.675 and 
zi = 0.39, respectively, and a flat FRW model is assumed. It should be noted that in 
reality there should be contribution from the eccentricity of the lens profile added to 
the above relation. Eq.(O-ll) and its counterpart for the NFW profile (Eq.(0-13)) in- 
dicate clearly that the lensing observable Be is primarily dependent on the lens model 
parameters, 9c and a, and less dependent on the cosmological parameters, Qm and 
w. (This sensivity was noted previously in Ref. for the case of an isothermal 
sphere.) This is essentially because strong lensing is sensitive only to the mass inside 
the Einstein radius. For example, in order to put a constraint on one needs to 
measure the Einstein radius and the velocity dispersion within 0{1)% accuracy and 
determine Qm within O(10)% accuracy. The former requirement would be an obser- 
vational challenge, while the latter could be accomplished. Another interesting 
observable is the location of radial arc, " which depends on the angular gradient 
of the projected mass, as expressed by Eq.(0-2). 

3. Constraining the lens model by measurements of the velocity dis- 
persion. We now show that there exists a degeneracy of lens models to a certain 
extent with regard to the projected mass density and that this degeneracy persists 
even when measurements of the velocity dispersion are included. We have in mind an 
arc or Einstein ring system, and thus one of the observablcs is the critical radius. Of 
course, the length of arc is another important observable, but we do not consider 
it here. 

As an illustration, we consider the well-known lensing system CL 0024-1-1654, 
although our argument is not limited to cluster lenses. Bright multiple arcs were 
discovered in the cluster CL 0024+1654 at z = 0.39 by Koo^^^ photographically. 
Five arcs are clearly seen in the HST image (see Fig. 1 in Ref. 23)). The redshift of 
the source galaxy was recently determined spectroscopically as Zg = 1.675.20) The 
distance to the arc from the center of the cluster is 9e = 34.6"^®), corresponding 
to 110/i~^kpc for the Einstein-de Sitter model. Here h is the Hubble parameter in 
units of lOOkm/s/Mpc. 
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Tyson et al. attempted to construct a high-resolution mass map of the cluster 
CL 0024+1654 using the Hubble Space Telescope. They found that the total mass 
profile within the arc radius is approximately represented by a power-law model, 



V-,/ \ _ -t*- V--- I 1-^ J 

^""^ ~ (l + x2)2-7' 



where x = r/vcore, K = TQOOilOO/iMopc-^, rcore = 35±3/t-^kpc and 7 = 0.57±0.02. 
They also noted that the asymmetry in the mass distribution inside the arcs for CL 
0024-1-1654 is very small (less than 3%). 

Recently, however, Broadhurst et al. have suggested that the mass profile 
of CL 0024-1-1654 is consistent with the NFW profile with Vg ~ mOh'^kpc and 
= Ps/Pcrit — 8000. Here Pcrit is the critical density. Similarly to Eq.(O-ll), for the 
case of the NFW profile, using Eqs.(0-5-0-6), we obtain the following relation for the 
above set of parameter values: 

^ c 2.0^ + 3.0^ - 0.23^ + 0.19^. (0-13) 

Shapiro and Iliev found that the projected mass density profile obtained by Tyson et 
al. is fit well by that obtained from the TIS profile with po — 0.064/i^Mqpc~^ and 
Tc — 20/i~^kpc. They also claimed that the mass profile obtained by Broadhurst et 
al. implies a velocity dispersion (> 2230km/s) that is much larger than the measured 
value. 

In Fig. 1, we show the projected mass density profiles for these three models. We 
assume the Einstein-de Sitter model and use two critical densities (pcritiz = 0) and 
Pcrit{z = zl)) for the fit with the NFW profile. We note that the angular resolution 
of the Hubble Space Telescope is 0.1", corresponding to 0.32/i~^kpc. The shaded 
region is the two-sigma interval of the mass profile determined by Tyson et al. We 
assume that the parameters {K,r core, 1) are Gaussian-distributed with dispersions 
equal to the error bars. Within the uncertainties on the fitting parameters, the 
power-law profile and the TIS profile look similar,*) while the NFW profile with 
Pcritiz = 0) deviates slightly from the power-law profile 

As is clear from Fig. 1, the problem with the NFW mass profile fitted by 
Broadhurst et al. may not result from a problem with the NFW model itself but a 
problem with the definition of pcrit used in the analysis. In the original NFW fit of 
the cold dark matter halo profile, pcrit in the relation ps = dcPcrit should be evaluated 
at the redshift of the object. However, within a fitting model, it is not necessary 
to do this, and we can treat ps just as a free parameter of the model. If we use 
Pcritiz = zl) as Shapiro and Iliev did, then the projected mass density is much 
higher than that obtained from the Tyson's fit of the data from the beginning. 

As is evident from the lens equation Eq.(OT), gravitational lensing provides 
information regarding the projected two-dimensional mass density. Therefore one 

Using the Davidon-Fletcher-Powell method, we independently fit Tyson's profile with the 

TIS model and found that Tyson's profile is also fit well (within lOOkpc) by the TIS profile with 
po ^ O-OSSTft^Mapc"^ and Tc ^ 15.4ft-^kpc. 

**^ In fact, T. Broadhurst informed us that they used pcrit{z = 0) to normalize ps. 
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Fig. 1. The projected mass density profile for CL 0024+1654. The dotted fine is the NFW profile 
with ps = 5cPcrit{zL) used by Shapiro and Iliev, and the short-dashed line is the same with 
ps ~ 5cpcrit{0), where, zl = 0.39 and we assume the Einstoin-dc Sitter model. The solid line 
the mass profile fitted by Tyson et al. with the 2a uncertainty by shades. The long-dashed line 
is the TIS profile fitted by Shapiro and Iliev. 



may wonder if the degeneracy of the lens model eould be lifted by combining with 
three-dimensional data, for example, the velocity dispersion. However, we believe 
that this is unlikely. The reason is the following. From Eq.(0-4) the sensitivity of 
the velocity dispersion on the parameters of the lens model is evaluated as 

— ~ + — . 0-14) 

Comparing with Eq.(O-ll), this indicates that the sensivity of the velocity dispersion 
is less dependent than the sensitivity of on the parameters of the lens model. The 
size of the system at which the velocity dispersion is measured is larger than the 
Einstein radius. 

In Fig. 2, we plot the circular velocity profiles divided by \/2 for each mass model, 
although it is known that v{r)/\/2 exactly coincides with the velocity dispersion 
only for a singular isothermal lens. The velocity profile of the power-law model was 
calculated using the Abel integral. As suggested by Fig. 1, the NFW fit with 
Pcritiz = Zl) predicts a velocity profile that is much higher than the measured value, 
in accordance with the claim of Shapiro and Iliev We find that the velocity 
profiles are all similar. The average velocity dispersion of CL 0024-1-1654 has been 
measured to be 1150km/s within a radius r ~ 600/i~^kpc, based on 107 galaxy 
redshifts , to an accuracy of roughly itlOOkms"-'^. For 33 galaxy redshifts, it has 
also been measured to be 1390km/s^^^. 
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Fig. 2. The circular velocity divided by ^/2, which coincides with the velocity dispersion for a 
singular isothermal lens, for the NFW model with ps = 5cPcrit{zL) (the dotted line), the same 
with ps = 5cpcrit{Q) (the short-dashed line) the power-law model (the solid line), and for the 
TIS profile (the long-dashed line). Here zl = 0.39 and we assume the Einstein-de Sitter model. 
The average velocity dispersion of CL 0024-1-1654 is measured to be 1150 km/s within a radius 
r ~ 600/i~^ km/s to an accuracy ± 100 km/s. 



4. Summary. We have examined the relation between lensing observables and 
model parameters (the lens model and the cosmological parameters) . We have found 
that observables are primarily dependent on the lens model parameters and have 
assessed the accuracy required to determine cosmological parameters [as expressed 
by Eq.(O-ll) and Eq.(0-13)]. It is the surface density of a lens that can be determined 
from the observations of gravitational lensing, and therefore there exists a degeneracy 
among lens models given the observational uncertainties. This degeneracy cannot 
be broken even by combining measurements of the velocity dispersion. We also 
have identified possible source of the problem of fitting the averaged mass profile of 
CL 0024-1-1654 with the NFW profile. The reconstruction of the mass profile using 
detailed shear maps with weak lensing observations may provide more accurate 

information regarding the mass profile of lenses. In any case, it is not possible to put 
meaningful constraints on cosmological parameters using gravitational lensing until 
we obtain more complete understanding of the lens model. 
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